% 假设数据服从均值为mu、标准差为sigma的正态分布
data = normrnd(mu, sigma, 1, 1000); % 生成1000个正态分布的随机数据
 
% 定义似然函数
logLikelihood = @(p, data) logpdf(Normal(p(1), p(2)), data);
 
% 使用mle函数进行最大似然估计
options = struct('Display', 'off'); % 关闭mle函数的显示
[estParams, ~, ~] = mle(data, logLikelihood, [0, 1], options);
 
% 输出估计的参数
disp(estParams);